A new correlation method to identify and separate charm and bottom production 

processes at RHIC 
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Electrons from semileptonic decays of heavy-flavor mesons (D and B) allow to study the energy 
loss of heavy-quarks in nuclear collisions at y/s = 200 GeV at RHIC. Since pQCD calculations have 
shown that the crossing point where bottom decay electrons start to dominate over charm decay 
electrons is largely unknown, an urgent need arises to access the relative contributions independently. 
A correlation method is proposed to identify and separate charm and bottom production processes 
on a statistical basis through tagging of their decay electrons and open charmed mesons. The 
feasibility for this method is demonstrated using PYTHIA and MCQNLO simulations. The latter 
allows to estimate the complete NLO contributions, including e.g. gluon-splitting diagrams. 

PACS numbers: 23.70.+j, 24.10.Lx, 25.75.Cj, 25.75.Gz 



INTRODUCTION 

Energy loss of partons is predicted to be a sensi- 
tive probe of the matter created in high energy nuclear 
collisions since its magnitude depends strongly on the 
color charge density of the matter traversed. In par- 
ticular, the understanding of the flavor dependent cou- 
pling of quarks and their fragmentation functions pro- 
vides key tests of parton energy-loss models and, thus, 
yields profound insight into the properties of the pro- 
duced highly-dense strongly interacting matter. Mea- 
surements at the Relativistic Heavy-Ion Collider (RHIC) 
at Brookhaven National Laboratory have revealed large 
medium-induced suppression at high transverse momen- 
tum (high pr) °f both the inclusive hadron yields and of 
back-to-back hadron pairs [1] . The principal energy loss 
mechanism underlying these effects is commonly thought 
to be medium-induced gluon Bremsstrahlung, which is 
expected to dominate collisional (elastic) energy loss for 
very energetic partons 0. 

Due to their large mass (to > 1 GeV/c 2 ), heavy quarks 
(charm and bottom) are believed to be primarily pro- 
duced by hard scattering processes (high momentum 
transfer) in the early stage of the collision and, therefore, 
are sensitive to the initial gluon density Q . Heavy-quark 
production by initial state gluon fusion also dominates 
in nuclear collisions where many, in part overlapping 
nucleon-nucleon collisions occur [4|]. Heavy-quark pro- 
duction by thermal processes later in the collision is low 
since the expected energy available for particle produc- 
tion in the medium (~0.5 GeV) is smaller than the energy 
needed to produce a heavy-quark pair (> 2.4 GeV). Theo- 
retical models based on perturbative Quantum Chromo- 
dynamics (pQCD) predicted that heavy quarks should 



experience a smaller amount of radiative energy loss in 
the medium than light quarks when propagating through 
the extremely dense medium due to the suppression of 
small angle gluon radiation 0, 0] . 

The energy loss of heavy-quark mesons is currently 
studied through the measurements of the pt spectra 
of their decay electrons. At high p?, this mechanism 
of electron production is dominant enough to reliably 
subtract other sources of electrons like conversions from 
photons and 7r° Dalitz decays. RHIC measurements in 
central Au+Au collisions have shown that the high pt 
yield of electrons from semileptonic charm and bottom 
decays is suppressed relative to properly scaled proton- 
proton collisions, usually quantified in the nuclear modi- 
fication factor (-Raa)@, S|- This factor exhibits an un- 
expectedly similar amount of suppression as observed 
for light-quark hadrons, suggesting substantial energy 
loss of heavy quarks in the produced medium. Energy- 
loss models incorporating contributions from charm and 
bottom do not explain the observed suppression suffi- 
ciently 0, [HI] • Although it has been realized that energy 
loss by elastic parton scattering causing collisional energy 
loss is probably of comparable importance to energy loss 
by gluon radiation [ll, 12 1, the quantitative description 
of the suppression is still not satisfying. Furthermore, 
it has been shown that collisional dissociation of heavy 
mesons in the medium may be significant in heavy-ion 
collisions [ll| ■ However theoretical models which include 
energy loss from charm only describe the observed sup- 
pression reasonably well [lfj . 

The observed discrepancy between data and model cal- 
culations could indicate that the B dominance over D 
mesons starts at higher px as expected. Theoretical cal- 
culations implying pQCD have shown that the crossing 
point where bottom decay electrons starts to dominate 
over charm decay electrons is largely unknown 
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Therefore, the relative contributions from charm and bot- 
tom meson decays to electrons have to be determined 
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FIG. 1: (Color online) Schematic view of the fragmentation 
of (a) a cc and (b) a bb pair. 



separately. 

This paper reports a new correlation method using az- 
imuthal angular correlations of heavy-quark decay elec- 
trons and open charmed mesons, which yields important 
information about the underlying production mechanism. 

CORRELATION METHOD 

In Quantum-Chromodynamics, flavor conservation im- 
plies that heavy quarks are produced in quark anti-quark 
pairs (cc and bb). A more detailed understanding of 
the underlying production process may be obtained from 
events in which both heavy-quark particles are detected. 
Due to momentum conservation, these heavy-quark pairs 
are correlated in relative azimuth (A</>) in the plane per- 
pendicular to the colliding beams, leading to the charac- 
teristic back-to-back oriented sprays of particles (dijet). 
A dijet signal appears in the azimuthal correlation distri- 
bution as two distinct back-to-back Gaussian-like peaks 
around Acf) = (near-side) and A<fi = tt (away-side). The 
correlation in their azimuthal opening angle survives the 
fragmentation process to a large extent in p + p colli- 
sions. Angular correlations of pairs of high p^ particles 
have successfully been used to study on a statistical basis 
the properties of the produced jets 

In this correlation method, charm and bottom pro- 
duction events are identified using the characteristic de- 
cay topology of their jets. Charm quarks predominantly 
hadronize directly to D° mesons (c — > D + X, BR = 
56.5 ± 3.2%) while bottom quarks produce D° via an in- 
termediate B meson (b -> B~ /W '/~B° -> D Q + X, BR 
= 59.6 ± 2.9%) [l6[. The branching ratio for charm and 
bottom quark decays into electrons is 9.6% and 10.86%, 
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FIG. 2: Transverse momentum spectrum of charm (solid line) 
and bottom decay electrons (dashed line) in (a) PYTHIA and 
(b) MC@NLO simulations of p+p collisions at yfs = 200 GeV. 



respectively. While triggering on the so-called leading 
electron (trigger side), the balancing heavy quark, iden- 
tified by the D° meson (D° -> A"~7r+, BR = 3.89%), is 
used to determine the underlying production mechanism 
(probe side). 

A charge-sign condition on the trigger electron and de- 
cay kaon provides a powerful tool to separate events with 
a cc or a bb pair. As an example, Figs. [T] (a) and (b) il- 
lustrate a schematic view of the fragmentation of a cc 
and a bb pair, respectively. Assuming the trigger lepton 
is an electron from the fragmentation of a c or b quark, 
the partner charm quark must be a c, hence producing 
a A~7r + pair. The bottom quark on the opposite side is 
a 5, which yield A + 7r~ pairs via the main decay mode 
B ->■ D 5 + X (BR = 59.6%). However, there is another 
channel, B -> D° + X (BR = 9.1%), which give R-TT+ 
pairs [16J. e~ 'K~ (e + K + ) pairs are also expected from 
semileptonic B decays, e.g., B~ — > D°e~17^. 

Thus, electron— kaon pairs with the opposite charge 
sign (called unlike-sign e — K pairs) identify B decays on 
the away-side of the azimuthal correlation distribution 
of decay electrons and D° mesons. Requiring like-sign 
e — K pairs select bottom on the near-side and charm 
and a small contribution from bottom (~15%) on the 
away-side of the e — D° correlation function. 

Requiring e — D° coincidence in the same event sig- 
nificantly improves the signal-to-background ratio over 
either technique individually. Moreover, the decay elec- 
trons provide an efficient trigger for heavy-quark produc- 
tion events. The shape of the azimuthal correlation dis- 
tribution allows a more differential comparison between 
the charm and bottom contributions owing to their dif- 
ferent decay kinematics. The feasibility for this correla- 
tion method is examined using PYTHIA and MC@NLO 
simulations. 
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PYTHIA V6.222 (a) 
B parent, trigger electron 





PYTHIA V6.222 (c) 
associated b— D° 




no trigger |) cut 

3<p <7GeV/c 
8<p T <20GeV/c 




FIG. 3: Transverse momentum distribution of the (a) B and 
(b) D mesons that yield trigger electrons in the indicated pr 
ranges. Panels (c) and (d) illustrate the transverse momen- 
tum distribution of the D° mesons from bottom and charm 
fragmentation, respectively, opposite the trigger electrons in 
the specified pr ranges. The electrons and D° decay products, 
kaon and pion, are selected having a pseudo-rapidity \r]\ < 1. 



MONTE CARLO SIMULATIONS 

The angular correlation function of charm and bottom 
decay electrons and D° mesons has been studied using 
leading-order PYTHIA simulations (version 6.222 with 
CTEQ5L PDF set, m c = 1.3 GeV/c 2 and m b = 4.5 
GeV/c 2 ) [HI of p + p collisions at ^ = 200 GeV. In 
total, 8 billion events are generated for charm (MSEL = 

4) and 4 billion events for bottom production (MSEL = 

5) with a cross section of 232 and 2.13 /xb, respectively. 
The default Peterson fragmentation function is used 
and the D/D* spin factor is taken into account [l8j]. 
Electrons within the pseudo-rapidity range of 1 77 1 < 1 
are assigned as trigger particles if they originate from 
charm (D° , D + , Df or their excited states) or bottom 
meson decays (B°, B + , B® or their excited states). The 
transverse momentum (pr) distribution of the decay 



electrons is shown in Fig. Ufa). Electrons from bottom 
decays starts to dominate over electrons from charm 
decays above pt > 4 GeV/c, consistent with results 
from pQCD calculations at the fixed-order plus next-to- 
leading log (FONLL) level [IlGjl. Figures E^a) and (b) 
depict the pt spectrum of B and D mesons, respectively, 
that yield trigger electrons in the indicated pr ranges. 
The associated D° mesons are accepted if their decay 
products (kaon and pion) fall within the pseudo-rapidity 
window 1 77 1 < 1. Figures [3jc) and (d) illustrate the pt 
distribution of the associated D° mesons from bottom 
and charm fragmentations, respectively. 

The azimuthal correlation function is calculated for 
all electron—!? and positron— D° pair combinations 
assuming a D° reconstruction efficiency of ^70% as 
typically observed in large acceptance experiments like 
the STAR detector [lj|]. In the following, we imply 
electron— D and positron— D° pairs when using e — D . 
Figures HJa) and (b) show the azimuthal correlation dis- 
tribution of heavy-quark decay electrons and D° mesons 
for like-sign e — K pairs from bottom production for two 
different trigger-electron px ranges. The same e — D 
correlation distribution is depicted in Fig. [4jc) and (d) 
for unlike-sign e — K pairs from bottom production and 
in Figs. 03a) and (b) for like and unlike-sign e — K pairs 
from charm production, respectively. Comparing the 
upper and lower panels of Fig. 0] one can conclude that 
like-sign e — K pairs select D° mesons from B decays 
on the near-side correlation whereas unlike-sign e — K 
pairs separate D° mesons from bb flavor creation on 
the away-side correlation. The near-side peak from B 
decays is relatively broad at intermediate px (3 < pt < 
7 GeV/c) and exhibits a double peak structure (cf. 
Fig. HJa)) which vanishes at higher px (cf. Fig.|4jb)). A 
comparison of the Figs. |D(a) and 03a) indicates that, for 
like-sign e — K pairs, the near-side peak is dominated 
by D° mesons from B decays whereas the away-side 
peak stems mainly from charm pair production (flavor 
creation) . The charm contribution for unlike-sign e — K 
pairs on the away-side is small (^14% compared to the 
like-sign e — K pairs) as shown in Figs. E)Jb). 



It has been shown [20} |2l| that higher order sub- 



processes like gluon splitting may have a significant con- 
tribution to the near-side correlation. The contribution 
from gluon splitting was determined using MCQNLO 
simulations of p+p collisions (version 3.3 with CTEQ6M 
PDF set) which allows modeling heavy-flavor hadro- 
production in a next-to-leading-order approach [22| . The 
MC@NLO computation uses the HERWIG event genera- 
tor (version 6.510) [23| for parton showering, hadroniza- 
tion and particle decays. 1 billion events are generated for 
each charm and bottom production with a cross section 
of 184 and 1.6 fib, respectively. The same particle se- 
lection criteria are used as for the PYTHIA simulations. 
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FIG. 4: Azimuthal angular correlation distribution of elec- 
trons and D° mesons from bottom decays generated in 
PYTHIA (solid line) and MC@NLO simulations (dashed line) 
of 200 GeV p + p collisions requiring like-sign (upper panels) 
and unlike-sign e — K pairs (lower panels) . The distributions 
are shown for trigger-electron transverse momentum ranges 
of (a+c) 3 < p T < 7 GeV/c and (b+d) 7 < p T < 20 GeV/c. 



Charm production 

like-sign e-K pairs 



-PYTHIA 
■ MCONLO 
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unlike-sign e-K pairs 




-PYTHIA J 




-MCONLO I 









FIG. 5: Azimuthal angular correlation distribution of elec- 
trons and D° mesons from charm decays obtained from 
PYTHIA (solid line) and MCQNLO simulations (dashed line) 
for (a) like- and (b) unlike-sign e — K pairs. Trigger-electron 
Pt range is 3 < pt < 7 GeV/c. Note the different scales for 
the correlation yield. 



agrees within 10-20% with the results from PYTHIA sim- 
ulations. This agreement is remarkable since these two 
event generators use different models for parton shower- 
ing and hadronization (fcj ordering in shower and string 
hadronization for PYTHIA and angular-ordered shower 
and cluster hadronization for HERWIG). The difference 
of the near-side peak in Fig. 03a) can be attributed to 
gluon splitting and is found to be (6.5±0.5)% of the open 
charm production observed in the studied pt range. 

Figure [6] depicts a two-dimensional plot showing the 
azimuthal correlation distribution of cc pairs (A<^(cc)) 
around the near-side peak of the azimuthal correlation 
distribution of e - D° pairs (A0(e, D )). The A0(cc) 
distribution exhibits a clear peak around zero which sup- 
ports the assumption that the near-side correlation peak 
of the A.<j>(e, D ) distribution is indeed from gluon split- 
ting. 



The pt spectrum of heavy-quark decay electrons is il- 
lustrated in Fig. IHb). Bottom decay electrons starts to 
dominate over charm decay electrons at a slightly lower 
p T compared to the PYTHIA results (cf. Fig.^a)). This 
seems to be due to the softer pt spectrum of the electrons 
from charm decays in the MC@NLO calculations. 

Figures H and [5] also show the results from MC@NLO 
simulations for the trigger normalized angular correla- 
tion function of electrons and D° mesons from bottom 
and charm production events, respectively. The corre- 
lation distribution from bottom production exhibits a 
similar shape as observed for PYTHIA simulations (cf. 
Figs. BJa-d)). The away-side peak shape of the corre- 
lation function from charm production (cf. Fig. Eta)) 



EXTRACTION OF THE RELATIVE BOTTOM 
CONTRIBUTION 

The relative bottom contribution for trigger electrons 
in the kinematical range 3 < pt < 7 GeV/c is obtained in 
two ways by comparison of the e — D° correlation yield on 
the near- (A(j) = ± ir/2) and away-side (A(j> = ir ± ir/2) 
from Figs. a+c) andEJa). 

Firstly by requiring like-sign e — K pairs which selects 
bottom on the near-side (cf. Fig. |4ja)) and charm on the 
away-side (cf. Fig. Eta)). The relative bottom contribu- 
tion — — is obtained from the D° yield on the near- 
side in Fig. [JJa) (D°(NS,b)) and away-side in Fig. Eta) 
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(D°(AS,c)) according to 



C-B 



1 



£)0(AS,c) 

_D°(NS,b) 



The branching ratio BR takes into account that D° from 
semileptonic bottom decays are always accompanied by 
an electron or more general by a lepton whereas electrons 
from charm decays have a probability of 56.5% to be 

ratio is found to 



balanced by a D° meson. The 
be 0.52±0.03 for PYTHIA and MC@NLO simulations. 

Secondly, the relative bottom contribution is deter- 
mined from the D° yield on the away-side which selects 
charm for like-sign e — K pairs (cf. Fig.[5ja)) and bottom 
for unlike-sign e — K pairs (cf. Fig.[4jc)). The c/b ratio 
is determined from the away-side D° correlation yield in 
Fig. Ha) (D°(LS, c)) and Fig. He) (£>°(ULS, b)) by 



£°(LS,c) 
L>°(ULS,b) 



c/b x 



BR(c^ D° + X) 
BR{b-> D°+X) 



PYTHIA and MC@NLO simulations give a, c/b ratio of 
1.01±0.07 and 1.27±0.09, respectively, for trigger elec- 
trons in the pt range 3 < px < 7 GeV/c. From 



('B 



1 



1 



c / b x §§F 



-*e+X) 



where the branching ratios for the c and b decays to 
electrons are quite similar, the e e + e is found to be 
0.53±0.05 and 0.47±0.04 for PYTHIA and MC@NLO, 
respectively. The uncertainties are obtained from the 
sum of the experimental uncertainties of the branching 
fractions in quadrature. 

The results obtained with the two different approaches 
agree within uncertainties. Furthermore, the extracted 
— — ratios show agreement with the relative bottom 
contribution from FONLL calculations 14, 1511 . 



SUMMARY 

The azimuthal angular correlation of heavy-flavor de- 
cay electrons and D° mesons in combination with a 
charge-sign requirement on electron and Z?°-decay kaon 
pairs allows, on a statistical basis, the separation of 
charm and bottom production and their sub-processes. 
The feasibility for this new correlation method is shown 
using PYTHIA and MC@NLO simulations which also 
yield an estimate of the complete NLO contributions (in- 
cluding gluon-splitting diagrams). The relative bottom 
contribution to the heavy-flavor decay electrons is deter- 
mined by comparison of the near- and away-side correla- 
tion distributions for charm and bottom production pro- 
cesses and is found to be ~50 % in the studied transverse 
momentum range 3 < pt < 7 GeV/c . 




A^(e,D u ) [rad] 

FIG. 6: (Color online) Azimuthal correlation of cc pairs 
around the near-side azimuthal correlation of e — D° pairs 
obtained from MC@NLO simulations. 
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